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Numerical simulations of the effect of a long-range scalar interaction (LRSI) acting only on non- 
baryonic dark matter, with strength comparable to gravity, show patterns of disruption of satellites 
that can agree with what is seen in the Milky Way. This includes the symmetric Sagittarius stellar 
stream. The exception presented here to the Kesden and Kamionkowski demonstration that an LRSI 
tends to produce distinctly asymmetric streams follows if the LRSI is strong enough to separate the 
stars from the dark matter before tidal disruption of the stellar component, and if stars dominate 
the mass in the luminous part of the satellite. It requires that the Sgr galaxy now contains little 
dark matter, which may be consistent with the Sgr stellar velocity dispersion, for in the simulation 
the dispersion at pericenter exceeds virial. We present other examples of simulations in which a 
strong LRSI produces satellites with large mass-to-light ratio, as in Draco, or free streams of stars, 
which might be compared to "orphan" streams. 
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I. INTRODUCTION 

Kesden and Kamionkowski ([7,^, hereafter KK06), in 
an independent application of the point made by Frieman 
and Gradwohl analyzed an important astronomical 
test of the weak equivalence principle, that the accelera- 
tion of a freely moving test particle is independent of its 
composition, applied to nonbaryonic dark matter (DM). 
They pointed out that if the dark matter significantly vi- 
olated the weak equivalence principle then the disruption 
of a satellite falling close to the host galaxy would tend 
to produce quite asymmetric stellar streams, a result of 
the non-inertial motion of the DM potential well that is 
holding the stars. This effect is contrary to the symmet- 
ric stellar streams of the Sagittarius (Sgr) satellite of the 
Milky Way (MW). The main purpose of this paper is 
to present a numerical simulation that demonstrates an 
exception to the Kesden-Kamionkowski effect. 

Our exception, symmetric stellar streams, requires two 
conditions. First, the departure from the equivalence 
principle is large enough that it separates the stars and 
dark matter in a satellite falling toward its host galaxy 
before significant tidal disruption of the luminous parts 
of the satellite. Second, the mass in the central region 
of the satellite is dominated by stars. This second con- 
dition is consistent with the observation that the stars 
in the Sgr galaxy and stream have heavy element abun- 
dance patterns [4] that are considered to be characteristic 
of a more massive satellite such as the Large Magellanic 
Cloud [5], where the mass in the luminous central parts 
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is dominated by stars [6]. We show in the simulation in 
jlTTl how the stellar system once separated from its dark 
matter halo suffers the usual pure gravitational tidal dis- 
ruption as it falls toward the denser parts of the host, and 
can end up with a symmetric stellar galaxy and stream 
with properties reasonably similar to the Sgr system. 

A significant departure from the weak equivalence prin- 
ciple in the dark sector would have other observable con- 
sequences. A simulation in §IVI is meant to approximate 
the evolution of Draco, which seems to be dominated by 
DM, perhaps because it is strongly enough bound and 
has remained far enough from the galaxy to have avoided 
complete separation of the stars and DM. Other simula- 
tions in this section show the development of a coreless 
stellar stream, which might be related to the Orphan 
streams in the MW [7], and twin gravitationally bound 
cores, one dark and one stellar. The latter has not been 
observed but might prove to be interesting. 

Our computation follows previous studies (reviewed in 
[8|] ) in expressing a departure from the weak equivalence 
principle as a fifth force of interaction in the dark sector. 
We take the interaction to be mediated by a scalar field 
that produces the fifth force, or LRSI, 

F = -^e-^/^=fi + lV, (1) 
r \r rj 

between DM particles. Here G is Newton's gravitational 
constant, m is the mass of a DM particle, r is the separa- 
tion vector, is a measure of the fifth force strength rel- 
ative to gravity, and is a "screening length" . This fifth 
force applies only between DM particles; it adds to the 
usual gravitational interaction among DM and baryons. 
In the simulations presented here is much larger than 
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the size of the galaxy so the fifth force is in effect an in- 
verse square law. In the scalar field model for Eq. ([1]) the 
screening length scales with redshift as ex {l-\-z)~^ 0. 
That is irrelevant for the effect on satellites of the MW 
but important in preventing any significant effect on the 
cosmological tests. 

KK06 showed that the ratio of stellar masses in leading 
and trailing streams of a disrupting satellite can be quite 
sensitive to the value of and they argued that the 
approximate symmetry of the Sgr stellar streams implies 
(3 < 0.04, in the notation of Eq. ([T|). Our alternative 
scenario for disruption follows if /3 is on the order of unity, 
meaning the DM fifth force is about as strong as gravity 
on the scale of galaxies. 

The next section outlines our numerical methods and 
parameter choices for building satellites of a MW-like 
galaxy and simulating the effects of gravity and an LRSI. 
The simulations use a version of the Gadget2 N-body 
code [10] modified to take account of the fifth force in 
Eq. ([T]). The numerical results for an Sgr-like satellite 
are presented in § IIIII Because the picture requires that 
the Sgr galaxy contains little DM we pay particular at- 
tention to the departure of the stellar velocity dispersion 
from virial equilibrium as the satellite moves around the 
galaxy. We argue in §|Vlthat this effect can reconcile the 
measured star velocity dispersion in the Sgr galaxy with 
a DM-free galaxy gravitationally bound by the mass in 
its stars. In § ITVl and IVl we comment on other issues to 
be explored in more detail. 



II. THE NUMERICAL MODELING 

We seek reasonably realistic simulations of the effect 
of an LRSI on observed satellites of the MW. In the sim- 
ulations the galaxy is represented as a static potential 
produced by a dark matter halo, a stellar disk, and a 
stellar bulge. The behavior of the satellites is modeled 
by the N-body dynamics of two particle species, DM and 
stars, with the stellar component initially concentrated 
well within the DM potential, as suggested by in 
initially spherical distributions. This description of ini- 
tial conditions and numerical methods is supplemented 
by details in Appendix [Al 



A. Mass Distribution in the Host Galaxy 

The gravitational potential of our model for the MW is 
represented by the sum of three fixed components. The 
first is a spherical DM halo with potential 



(2) 



where r is the distance from the center of the galaxy, 
'^haio = 131.5 km s~^, and b = 12kpc. The second is a 



Miyamoto-Nagai disk [l2| with potential 

GMdisk 



disk 



(3) 



where z and r± are cylindrical coordinates and Mdisk = 
10^1 Mo, a = 6.5 kpc, and d = 0.26 kpc. The third is a 
spherical Hernquist bulge [isj with potential 



GMbulge 



(4) 



where Mbuige = 3.4 x 10^° M© and c = 0.7 kpc. 

The DM particles feel the additional force given by the 
gradient of <I>scaiar = P^haio- We neglect the exponential 
term in the scalar force in Eq. ([T]) because the size of the 
galaxy is supposed to be much smaller than the screening 
length r^. 



B. Construction of the satellites 

The initial DM and star distributions in a satellite are 
constructed to be spherically symmetric. The virial ra- 
dius is r^.^ = r2oo/A, where r^^^ is the distance within 
which the mean mass density of DM plus stars is 200 
times the present critical cosmological density (for Hub- 
ble parameter Hq =70 km s~^ Mpc~^), and A is dis- 
cussed in the next subsection. The DM density profile at 
^ ^ ^vir Navarro, Frenk & White form [l3, flBj, 



Pdm(0 
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r/r^.^ (1 + c^r/r^.J 
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where is the concentration parameter and po is a nor- 
malization factor. At r > r ^.^ we choose an exponential 
DM density profile, as in [16|], 



exp 



The continuity of dp^^/dr at r = r^.^ demands 



-1 - 3c 
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The undisturbed inner stellar component follows the 
modified Hubble density profile [17], 



Pb(0 =Pi 



-2/3 



at r < , 



(8) 



where is a core radius, pi is a normalization factor, 
and is a truncation radius of the stellar component. 
At r > the stellar density run is 



-3/2 



— exp 



(9) 
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where 

«b =rj- , (10) 

from continuity of dp^/dr at r = r^. 

We assign the satehite particles — DM and stars — 
statisticahy isotropic velocity distributions that produce 
near dynamical equilibrium prior to significant perturba- 
tion by the host galaxy, using the variant of the method in 
[l^ outlined in Appendix [Al Before placing the satellite 
in the gravitational potential of the host galaxy we com- 
pute the motions of the satellite particles for 20 gravita- 
tional dynamical times (as measured at the scale radius 
r^.Jc^) using the Gadget2 N-body code. The evolved 
density profiles are quite stable over this range of time. 



C. Parameter choices 

Throughout this study we set the ratio of the total 
DM to stellar mass within a satellite to be close to the 
global cosmic value [l^ . The stellar mass fraction may be 
smaller than that in a satellite, as it is in larger galaxies 
(as reviewed in [17]), but that need not matter because 
excess DM in the outermost parts of the satellite could be 
stripped away before much happens to the stellar com- 
ponent. 

We summarize in table |I] the parameters for the initial 
structures of the six model satellites we have chosen to 
illustrate possibly interesting effects of an LRSI. Entries 
1 to 6 are the scale lengths in kpc and 7 and 8 are the 
concentration parameters defining the structure of the 
model satellites, as defined by Eqs. ([5]) to ([10]), prior to 
disturbance by the host. 

The stellar component is embedded in a generally more 
extended DM component. The parameter = r^/r^^ 
in entry 9 is the ratio of the stellar to DM core radii. 
According to [11], the range of values of this parameter 
in MW satellites is quite large, ranging from perhaps 
~ 0.05 to close to unity. The strength of the LRSI 
relative to gravity is listed in entry 10. The simulation 
SI includes gravity alone (/3 = 0). 

The host galaxy is centered on the origin of the coor- 
dinates, the disk is centered on the Z = plane, and the 
satellite moves in the YZ plane. The initial satellite po- 
sition is in the plane of the disk at the distance (in kpc), 
from the center of the galaxy, given in entry 11. The 
initial velocity in entry 12 (in km s~^) is in the positive 
Z direction. 

Entries 13 to 22 give the initial stellar and DM masses, 
in units of 10^ M©, within various of the radii listed in en- 
tries 1 to 4. The overdensity used to set the virial radius 
is 200 X A^, where A is listed in entry 23. The exception 
to A = 1 is S3, where A = 2.7. This satellite is sup- 
posed to correspond to a high concentration (c^ = 13.5) 
DM halo. However, the number of particles needed for 
proper numerical modeling is proportional to the fourth 
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TABLE L Simulation parameters. See description in text. 

power of the concentration (see Appendix [A|) . This be- 
comes numerically expensive, and, as a compromise, we 
work with a larger A to imitate the desired larger con- 
centration. 

The parameters N^^^ and N^^ are the numbers of parti- 
cles, in units of 10^, representing the mass distributions 
in stars and DM. The parameters e^^^ and e^^ are the 
lengths in parsecs of the smoothing of the inverse square 
law for DM and stars in the numerical simulations, as 
defined in the Gadget2 code. 

Six simulations for the initial conditions in TableUhave 
been run with Gadget2. Only SI does not include LRSI; 
it serves as a benchmark for comparison to the standard 
picture. Simulation S2 has initial conditions similar to SI 
but is run with LRSI at the strength of gravity (/3 = 1). 
It is intended to produce the general properties of the Sgr 
galaxy and stream. The satellite S3 is constructed with 
5^ = 1, and has a DM-dominated core. It becomes almost 
completely dominated by DM as a result of an LRSI. It 
could be an analogue of Draco, though it is not intended 
to be a close match to any of its properties other than 
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FIG. 1: The initial stellar and DM masses as functions of 
distance r from the center of the satellites SI and S2. 



the high M/ L and the apparent absence of tidal streams. 
In S4 the stellar core is destroyed while the DM core 
survives, leaving a smooth stellar stream. Simulations 
S5 and S6 are attempts to illustrate the formation of 
twin cores, DM and stellar. They share the same initial 
conditions but differ by the LRSI strength jS. For the 
Gadget2 code parameters, we chose conservative values, 
resulting in conservation of at least 99.5% of the total 
energy for all simulations. 



III. MODELING THE SAGITTARIUS SYSTEM 

We describe here the evolution of model satellites SI 
and S2 that are meant to simulate the formation of the 
Sgr galaxy and stream without and with the effect of 
a strong departure from the weak equivalence principle 
represented by an LRSI in the dark sector. We compare 
the models to observations in ^ 

Initially the cores of SI and S2 are dominated by the 
stellar mass out to 4 kpc radius, as illustrated by the DM 
and stellar mass density runs in Fig.[Tl As we have noted, 
and will discuss in more detail in ^ this is consistent 
with the fact that the Sgr stars have been compared to 
the population in the Large Magellanic Cloud (the LMC; 
e.g. [5]), in which the evidence is that the mass in the 
central region is dominated by stars 0. 

To get similar orbits we have assigned the same initial 
velocity but slightly different initial positions for SI and 
S2 in the host galaxy (entries 11 and 12 in table [Tj). Both 
simulations are run for 6 Gyr. 

Evolutions of the particle distributions projected in the 
Y Z plane perpendicular to the plane of the galaxy are 
shown in Figs. [2] and [3l The stellar core in SI remains 
embedded in its DM halo as the satellite develops quite 
symmetric tidal streams. Early on the streams contain 
mainly DM, but after the first pericenter passage stars 



become prominent in the streams. The evolution of S2 
is quite different, but it ends up with similar symmetric 
stellar streams. In a frame of reference moving with the 
DM core, the stellar core experiences an effective force 
that works against the gravitational pull of the DM core. 
At = 1 and at the initial position of S2 this effective 
force overcomes the gravitational pull of the DM core. 
The result is a nearly complete segregation of the stars 
and DM well before the first pericentric passage. The 
separated stellar core is then free to develop symmetric 
stellar tidal streams under purely gravitational dynamics. 
Meanwhile the DM core is completely disrupted. 

The behavior of S2 is very different from the KK06 
simulations. They focus on cases where the strength of 
the LRSI is weak enough that there is not segregation of 
the stellar and DM components prior to significant tidal 
disruption of the stellar system. In the KK06 simulations 
the effective force on the stars in the DM halo rest frame 
tends to pull stars out of one side of the halo, producing 
quite asymmetric stellar streams. 

Fig. m shows measures of the evolution of SI and S2. 
The vertical scales of each curve are arbitrarily normal- 
ized to fit the figure, independently in each panel. The 
dot-dashed curves are the satellite distance from the cen- 
ter of the host, showing pericenter passages. The solid 
blue curves show the rate of loss of the satellite stel- 
lar mass, and the dotted curves show the residual stellar 
mass. Here, the satellite is defined as the set of self-bound 
particles, which for SI is calculated with both DM and 
stellar particles, and for S2 with stellar particles alone. 
As expected, the mass loss rate in both simulations is 
highest during pericenter passages where the tidal force is 
strongest. The figure also shows one-dimensional stellar 
velocity dispersions measured within 0.7 kpc projected 
distance from the center of the satellite. The line-of-sight 
velocity dispersion for an observer at (0, — 8.5kpc, 0) rel- 
ative to the center of the host galaxy is plotted as the 
thick black lines. The red dashed lines show the disper- 
sions in the directions giving the maximum and minimum 
values. In these simulations the line-of-sight direction is 
close to the direction of maximum velocity dispersion at 
pericenter passages. 

The satellite is constructed to be in dynamical equi- 
librium with an isotropic velocity distribution prior to 
disturbance by the host. The anisotropy that develops 
during pericenter passages, and the increase in the line- 
of-sight velocity dispersion, are combined results of the 
contamination of escaping particles and the transient de- 
parture from dynamic equilibrium, including the distor- 
tion of the satellite shape caused by the strong tidal force. 
The deviations from an isotropic stellar velocity distribu- 
tion become larger in S2 once the stars are separated from 
the DM, the residual DM in SI keeping the stars more 
tightly bound. 

Fig. [5] illustrates how transient disturbances to the stel- 
lar velocity distribution can affect dynamical mass esti- 
mates of the satellites SI and S2. The blue curves show 
the evolution of the stellar mass and the stellar line-of- 
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FIG. 2: Snapshots in the YZ plane of the evolution of the SI satellite (without an LRSI). The stellar and DM particles are 
plotted as red open circles and blue dots, respectively. For clarity, only a subset of the stellar particles are plotted. The yellow 
star indicates the location of the most bound stellar particle. 




FIG. 3: The same as Fig. [2] but for S2 with an LRSI with p = 1. 
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FIG. 4: Measures of evolution of the SI (panel to the right) and S2 (to the left) simulations. 
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FIG. 5: Evolution of the satellite bound stellar mass and stellar line-of-sight velocity dispersion follows the blue lines commencing 
near the top right corner of each panel and ending toward the lower left. At the thick portions the velocity anisotropy is less 
than 15%. The dashed line varies as cr^^, as expected for dynamical equilibrium in a satellite with fixed size. 



sight velocity dispersion a^, both computed within pro- 
jected radius 0.7 kpc perpendicular to the line of sight. 
The thick portions of the curves correspond to times 
when cr^^^^ - a^^^^ < O.lSa^^^^^, where a^^^^ is the 
RMS value averaged over all directions. The satellites 
start at the top-right corners and end toward the lower- 
left after 4.17 Gyr for SI and 3.3 Gyr for S2. The dashed 
red lines show how the velocity dispersion would vary 
with the mass, M(< r) ex a^^, at dynamical equilib- 
rium and fixed radius r. The blue curves pass close to 
this condition when the velocity anisotropy is small. At 
the thin portions of the line the satellite is closest to 
the host galaxy and most disturbed. Here the mass is 
smaller than expected from the velocity dispersion 
under the assumption of dynamical equilibrium, by a 
factor as large as two and even larger for S2 towards 
the final time. Earlier studies of the possible importance 



of tidally-driven departures from dynamical equilibrium 
differ, but some indicate that velocity dispersions of long- 
lived dwarf satellites can be interpreted as the effect of a 
departure from equilibrium rather than the presence of 
dark matter [2(51, 12J, [H, [13 . 

Fig. [6] shows that the stellar velocity dispersions in the 
cores of SI and S2 are close to constant across the face 
of the galaxy. The curves at t = 3 Gyr are higher than 
at t = 2.8 Gyr in both SI and S2, because the former 
coincides with a pericentric passage in both simulations. 

The evolution of the mean surface density profile as 
a function of projected distance R from the center of 
the satellite is shown in Fig. [7l The stellar component 
plotted as the thin lines preserves a flat density in the 
inner regions, in agreement with [IJ] 

The parameter KK06 use as a measure of the degree of 
symmetry of the tidal stream is the ratio of the number of 



7 



20 
18 

" 16 









\ /X 


- -X- - - ^ P=0, t=3 Gyr 








. x X P-0, t-2.8 Gyr 




•^""^ P=l,t=3Gyr 




^ X P=l,t=2.8Gyr 


10-^ 


10^ 




R [kpc] 



FIG. 6: Line-of-sight stellar velocity dispersions as functions 
of projected distance from the satellite center are shown for 
SI (/3 = 0) and S2 (/3 = 1) at various times indicated in the 
figure. 



leading to trailing stars. KK06 find that in all their simu- 
lations without an LRSI this ratio is never lower than 0.5, 
while in the simulations with 0.04 < /3 < 0.09 it never 
exceeds 0.2. Indeed, this supports the conclusion that a 
weak LRSI is seriously constrained by the observations of 
Sgr. In the strong LRSI case we are considering, there is 
some asymmetry in the early stages of evolution, which 
begins to decay right after the DM leaves the stellar core. 
At t = 1.7 Gyr the ratio of the number of leading to trail- 
ing stars, excluding particles that are closer than 4 kpc to 
the center, is 1.09 and 0.817 for SI and S2 respectively. 
At t = 2.85 Gyr the ratios are 0.96 and 0.88. 

We discuss how these models for the Sgr galaxy and 
stream compare to the observations in ^ after present- 
ing simulations of a few other situations of possible as- 
tronomical interest. 



IV. OTHER MANIFESTATIONS OF AN LRSI 

The mass-to-light ratio of the Draco satellite, at dis- 
tance D ^ SO kpc from the MW, is estimated to be 
M/L > 100 Solar units [25]. This large value is not 
reasonably explained by deviations from equilibrium [2^ 
such as those occurring near pericenter passages dis- 
cussed in the previous section; it more likely means the 
mass is dominated by DM throughout this satellite. In 
the standard scenario, tidal mass loss may result in an 
increased M/ L ratio [24] , however this effect is not strong 
enough to explain highly DM dominated satellites as 
Draco. LRSI offers a mechanism that can work either 
separately or in tandem with the usual postulate that 
photoionization heating by the UV background has sup- 
pressed star formation [2^, [2^ , provided the LRSI is 
strong enough to have stripped away most of the stars 



but weak enough to leave the most tightly bound stars 
in a DM halo that remains bound. 

Our simulation S3, with /3 = 1, illustrates a case of 
significant but not complete loss of stars in the LRSI pic- 
ture. The stellar and DM core radii initially are the same 
{s^ = 1 in tablelj), and the central region is dominated by 
DM (at close to the cosmic baryonic to DM mass ratio). 
One sees in Fig. [8] that the LRSI pulls most of the stellar 
particles out of the DM halo early in the evolution, in the 
direction opposite to the motion. That leaves a bound 
DM halo with a tighter concentration of stars, in which 
the ratio of bound DM mass to stellar mass is 16 times 
larger than the initial value. 

Simulation S4 illustrates the development of the core- 
less stellar stream in Fig. [9] that might be compared to 
the orphan streams in the MW. In this simulation the 
stellar core is not gravitationally self-bound (the satellite 
core mass is not dominated by stars) so it is disrupted 
by the segregation driven by the strong LRSI. Even af- 
ter 7.4 Gyr the stars remain in a smooth stream, little 
disturbed by the small mass left in the DM core in this 
simulation. 

Simulations S5 and S6 are tuned to explore the possi- 
bility that both the DM and stellar cores survive sepa- 
ration by an LRSI, ending up as separate self-bound ob- 
jects. Since a DM core is bound by the combined effect of 
gravity and the LRSI its binding energy is enhanced by 
the factor l-\-/3 relative to a stellar system with the same 
mass distribution. This means that, once we arrange for 
a self-bound remnant stellar core, a DM core will survive 
if /3 is large enough. We did not find bound twin cores 
in S5 with (3 = 1. Fig. [TOl for /3 = 2 (S6), shows twin 
cores, but the remnant bound DM halo mass is only 1% 
the mass of the original satellite. 



V. DISCUSSION 

The central point of this paper is that a strong depar- 
ture from the weak equivalence principle for nonbaryonic 
dark matter, with P c:^ 1 in the LRSI model in Eq. ([1]), 
can produce a system that resembles the Sgr galaxy and 
stream from a progenitor that seems reasonable from the 
point of view of the astronomy. We have also taken note 
of other modes of satellite disruption that might be ob- 
servationally interesting. 

It is important for our picture of the origin of the Sgr 
system that the mass distribution in the central part of 
the Sgr progenitor is dominated by stars. The thought 
that this might be so is motivated by the observation 
that stellar element abundances in the Sgr galaxy and 
stream have been compared to LMC stars (e.g. [5]), and 
the evidence that the mass in the luminous central parts 
of the LMC is dominated by stars Q. 

We have aimed for a progenitor in our simulation that 
has structure similar to the LMC (though a rough sim- 
ilarity is likely all that could be meaningful). Our SI 
and S2 cases each have initial total mass ~ 2.2 x 10^ M©. 



8 




10 
R [kpc] 



FIG. 7: The mean stellar surface density profiles at different times for SI (right panel) and S2 (to the left). The DM profiles 
for SI are plotted as the thicker lines. 
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FIG. 8: Illustration S3 of the development of a satellite that retains a small but nonzero fraction of its stars. The orientation 
is the same as in Figs. [2]and[3l Red open circles are the stellar component, blue dots the DM. 



The virial radius is slightly less than 20 kpc, and stars 
dominate the mass within 4 kpc radius. Kim et al. ^] 
find that the LMC mass is ^ 3.5 x 10^ M© within 4 kpc 
radius, and the total disk mass is 2.5 x 10^ within 
7.3 kpc radius. By these measures the initial structures 
of SI and S2 are reasonably similar to the LMC. 

In S2, with the strong LRSI, the bound mass in stars 



at the third pericenter passage at 3 Gyr is 5.2 x 10^ Mq. 
At stellar mass-to-light ratio M/L = 2 Solar units the 
luminosity would be L ~ 1 x 10^ Lq. This is fairly close 
to the measured luminosity of the Sgr galaxy, in the range 
(2-5.8) X lO^L© 13, 

In our model the Sgr galaxy no longer contains dark 
matter. Estimates of M/L for this galaxy are larger than 
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FIG. 9: Mass distributions near the end of simulation S4. 




FIG. 10: Simulation S6. Circled in red is the remnant DM core. 



expected for stars; thus from a recent program of obser- 
vations [32[ finds M/L = 17.5. This assumes dynamical 
equilibrium, however, and we have illustrated in Figs. [H 
and [5] how the galaxy may be disturbed from equilibrium 
by its host. In S2 the satellite at 3 Gyr is at pericenter 
18 kpc from the center of the host galaxy, similar to the 
present distance of the Sgr galaxy, 16 kpc [33]. Also, al- 
though our host galaxy is rigid its mass distribution is 
designed to resemble that of the MW. At this third peri- 
center passage the line-of-sight velocity dispersion within 
the innermost 0.7 kpc of the simulated satellite is close to 
constant across the core at cr = 9.8kms~^ (Fig. [6]). This 
is in line with the finding of [32[ that the velocity disper- 
sion in the Sgr galaxy is close to constant across the inner 
0.7 kpc from the center defined by M54 at a = 9.6 kms~^. 

The departure from dynamical equilibrium at pericen- 
ter, and the large apparent value of the mass, also occurs 
in conventional gravity without the LRSI, as in SI, and 
in the simulations by [sl, [35[. The effect is larger in S2 



than SI, however, because of the greater fragility of the 
stellar core after the early loss of the DM. 

Our conclusion from these comparisons of S2 to the 
observations is that it is possible to obtain a reasonably 
realistic picture of the origin of the Sgr system, including 
a DM-free galaxy and a symmetric stream, under the 
effect of a strong departure from the weak equivalence 
principle in the dark matter. It would be interesting to 
study the effects of LRSI on the features of the tidal 
streams, as done in [36]. 

Our simulations suggest alternative interpretations of 
other observations. The example in Fig. [HI a satellite that 
loses most but not all of its stars, might be compared to 
Draco. A rough analytic argument that the Draco DM 
halo can retain its most tightly bound stars against the 
effect of a strong LRSI goes follows. If the MW rotation 
curve for baryons is close to flat at Vc = 220 km s~^ at 
distance D then a star in the DM halo rest frame expe- 
riences acceleration f3v'^/D. The gravitational accelera- 
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tion produced by the mass in the satehite at distance r 
from its center is maximum at GM(< r)/r'^ :^ 2cr^/r^ 
at the core radius r^. According to [25[ the stehar core 
radius of Draco is 160 pc. These numbers indicate that 
the one-dimensional stehar velocity dispersion must be 
a = 7/3^/^ km s~^ to keep stars in the core of Draco at 
its present distance D. Since the measured dispersion is 
in the range a = 8.5 — 10.7 kms~^ [25], we conclude that 
a strong LRSI, ~ 1, allows Draco to contain stars as 
long as it has not ventured much closer to the MW than 
its present distance. 

In the numerical simulation in Fig.[8]there are remnant 
bound stars, the ratio of bound DM to stellar mass is an 
order of magnitude larger than the original value with the 
DM remaining dominant at all radii, and the stars that 
have escaped are broadly scattered. This is in line with 
the properties of Draco. Whether more detailed measures 
of the simulated remnant satellite can be consistent with 
what is observed in Draco is a subject for separate study. 

Another example for closer study is the difference be- 
tween the observed numbers of satellites and the far 
greater numbers of DM halos predicted by the ACDM 
cosmology [s^, S, |39|. This usually is taken to mean 
that many DM halos are not visible because they were un- 
able to accrete or retain photoionization-heated plasma 
[27|, [28|, |29|. The possibility offered by an LRSI is that 
many halos are not visible because stars have been re- 
moved from their DM halo, perhaps largely, as may have 
happened to Draco, or entirely. 

Depending on the structure of the progenitor, stars re- 
moved from the DM may be broadly dispersed or may 
end up in a coreless stream, as in Fig. [9l The latter 
offers an interpretation the orphan streams in the MW 
(e.g. Orphan streams can form by tidal disruption 

in standard gravity acting on satellites in pericentric pas- 
sages close enough or repeated often enough to make it 
difficult to identify the remnant satellite or totally dis- 
rupt it. In this scenario the formation time of coreless 
streams could be quite long [g^. In the LRSI scenario 
demonstrated in Fig. [9] coreless streams may form rapidly, 
and could be found in orbits with large pericentric dis- 
tances. This behavior might be detectable, and is a sig- 
nature of a strong LRSI. 

Some stars separated from the DM may end up in grav- 
itationally bound DM- free dwarf galaxies. There are dis- 
cussions of dwarfs with mass-to-light ratios that might be 
considered characteristic of a normal stellar population 
(4Q| . DM- free dwarfs would be an interesting challenge to 
the standard theory. If the stars in the progenitor were 
much more strongly concentrated than the DM then a 
normal gravitational tide could strip away most of the 
DM, leaving a bound system dominated by stars, but it 
would be a delicate operation. An LRSI offers an alter- 
native that operates in a straightforward way. 

A recent study [4l[ challenges the model presented 
here. The main concern is that the satellite S2 is tightly 
constrained in the initial distance in which it may first en- 
ter the MW, since it is too light for dynamical friction to 



modify its orbit into the observable one. It is important 
to notice that in the standard model (without LRSI), Sgr 
could have many different possible orbital histories which 
lead to its currently observed setup. According to [42^, in 
one family of solutions the Sgr may have entered the MW 
at a galactocentric distance close to its current apocentre 
of 60kpc, with a mass of 10^ M©. In this model, the Sgr 
is indeed too light for gravitational friction to have any 
significant effect on its orbit. This model is consistent 
with our satellite model S2. In another distinct family of 
solutions, the Sgr could enter the MW at a distance of 
more than 200 kpc with a total mass of about 10^^ M© 
which is 100 times more massive than S2. In this sce- 
nario, dynamical friction plays a key role in bringing Sgr 
to its current orbit, while along the way it looses more 
than 99% of its mass. This scenario is consistent with 

0. 

We have run a simulation of this 'heavy' Sgr scenario, 
including a 'live' MW halo to account for dynamical fric- 
tion. The results are illustrated in Fig.[TTl Here, the MW 
halo is modeled with an NFW profile of concentration 
15, mass of 2 x 10^^ M©, and virial radius 250 kpc. The 
number of particles is 2 x 10^, and the smoothing 8 kpc. 
The Galactic disk and bulge have the same parameters 
as in the previous simulations, however the bulge is not 
a static potential, but a single massive particle smoothed 
at 9 kpc with a spline kernel, as all particles in Gadget2, 
and the center of the stellar disk static potential follows 
the particle that represents the bulge. The initial Sgr 
has a baryonic component of mass 4 x 10^^ M© with a 
truncation radius of 3.6 kpc and is represented by 7 x 10^ 
particles, smoothed at 55 pc. The DM component is five 
times more massive, has a virial radius of 29 kpc repre- 
sented with 2 X 10^ particles, and smoothed at 850 pc. 
Initially, the Sgr is positioned in the stellar-disk plane 
at a GC distance equal to the Galactic virial radius of 
250 kpc, and its velocity is set to the baryonic circular 
velocity of 185kms~^, perpendicular to the stellar-disk 
plane. 

At the initial position in this simulation the mass and 
concentration of the Sgr DM halo are large enough to 
prevent separation between stars and DM. This is im- 
portant, because while the baryons are bound to the DM 
core, the orbital dynamics is almost independent of (3, 
since changing /3 is equivalent to a change in the gravita- 
tional constant G, which in turn is equivalent to a trans- 
formation in the time variable t ^ t x ^/^~^^. Hence, 
the Sgr begins to sink in the usual way, while experienc- 
ing dynamical friction and loosing mass to the MW. Only 
when the Sgr is in an orbit similar to the one which is 
inferred directly from the observations (i.e. with an apoc- 
enter of 60 kpc), the DM halo becomes small enough to 
allow full separation. Afterwards the Sgr develops sym- 
metric tidal arms, in a similar way to SI and S2, and 
because it looses about 99% of its mass, it is similar to 
the Sgr today. This works also for larger turn-around 
distances, since as mentioned above, the dynamics are 
independent of P as long as the DM and stellar cores are 
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initially tightly bound. 



VI. CONCLUDING REMARKS 

We have addressed some of the issues in astronomy 
that are raised by the postulate of a strong departure 
from the weak equivalence principle by nonbaryonic dark 
matter. One must also consider the effect on structure 
formation on larger scales. A preliminary analysis [44] 
of the effect of adding a strong LRSI to the ACDM cos- 
mology indicates an improved fit of theory to observa- 
tions of voids, which are more completely emptied by 
the LRSI, and of giant galaxies, which complete major 
mergers sooner. Since the DM clusters more rapidly in 
the LRSI scenario, the baryons on scales below r^^ would 
be expected to lag behind the DM. This effect could leave 
a fraction of the baryons out of galactic haloes and help 
explain the missing baryon problem reviewed in [45[ . We 
plan to present a more detailed and systematic analy- 
sis of these considerations using the numerical methods 
developed here with the further adjustment of the nu- 
merical code to take account of the departure from an 
LRSI inverse square law in Eq. ([1]). 

As we noted in the last section there are still other 
issues to consider. One is whether a satellite with all 
the properties of Draco may naturally form in the strong 
LRSI picture. Another is the excess satellite problem 
[37I . [ssL [3^ . We stress here that LRSIs do not necessarily 
imply DM free satellites. In fact, small but high concen- 
trated DM halos could easily trap the stellar components 
so that full segragation never happens. This means that 
LRSI may even help increase the M/L ratio by yanking 



only the least bound stars. An interesting quantitative 
test in a pure particle simulation of the formation of an 
galaxy would count satellites that are entirely stripped 
of the particles that represent baryons, satellites that re- 
tain enough stellar particles to be visible, and satellites 
that fall in at low redshift and are capable of producing 
symmetric and orphan stellar streams. 

We are considering a substantial departure from stan- 
dard physics, but it is in the little explored dark sector 
where our present physics is so very simple as to seem 
possibly suspicious. Thus we feel that explorations of 
more complicated dark sector physics, while speculative, 
are worthwhile. Our studies of the observational tests of 
the speculative idea that there is a long-range departure 
from the weak equivalence principle in the dark sector 
have not proved to be discouraging so far. 
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APPENDIX A: SATELLITE CONSTRUCTION 

We follow the recipe in [l^, with minor adjustment to 
accommodate two components, stars embedded in a more 
massive DM halo. In this method, the single particle 
distribution function / in phase space is a function only 
of the dynamical constant Q, 



2rJ 



(Al) 



where e is the total binding energy per unit mass, L is the 
particle angular momentum vector per unit mass, and 
is the scale of anisotropy, given as 



1 



(Jr + Va 

The distribution function gives the number density. 



(A2) 



P{r)= lf{s,L)d'v. (A3) 



The inverse of this equation is [4£ 



fiQ) 



f d^PQ J_ f dpQ 
J d^^ ^ ^\d^ 



(A4) 



where pq = p{r){l + r^/ro^) and = — $ — is the 
relative potential, including its gravitational and scalar 
components. 

We construct two distribution functions for the DM 
and the stars in a satellite. The relative potential of the 
stars is 



G 

-[M,{r) + M,^{r)]dr, (A5) 



where (r) and M^^ (r) are the masses in stars and DM 
within radius r. The relative potential of the DM is 



dm 



G 

^ [M, (r) + (1 + /3)M,„ (r)] dr . (A6) 



Convergence criteria [47|] for numerical sampling of the 
mass and velocity distribution functions require that the 
number of particles within the virial radius satisfies 



N. > 



(ln(l + cj-cj(l + cj)2 



(A7) 



for an NFW profile with concentration c^, and the 
smoothing length e is 



6 = 4- 



(A8) 
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Under these conditions there are enough particles that 
the relaxation time is much longer than the simulation 
dynamical time, and particle accelerations do not exceed 
the limit set by the smoothing length that keeps the 
acceleration induced by two-particle interactions weaker 
than the mean field halo acceleration. 

Our numerical satellite construction code applied to 
a Hernquist potential [13] yields a distribution function 



consistent with the analytic solution to better than a part 
in 10^ except at the edges of the distribution. We find 
similar consistency for the combined NFW and plummer 
profiles, where we numerically constructed the distribu- 
tion functions from the analytic profiles and then numer- 
ically reconstructed from the distribution functions the 
density runs of the two components. 



